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Sodalite,  N.  P.  Blake  and  H.  Metiu,  J.  Chem.  Phys. ,  submitted 

2.  Electron  Traps  and  Polaron  Formation  in  Dehydrated  Sodalite,  N.  P.  Blake,  P.  C.  Weakliem, 
and  H.  Metiu,  J.  Chem.  Phys.,  submitted 

3.  knAb-Initio-Qd&Q&  Transferable  Potential  for  Sodalites,  N.  P.  Blake,  P.  C.  Weakliem,  and  H. 
Metiu,  J.  Phys.  Chem.,  to  appear 

4.  Fitting  Potential  Energy  Surfaces:  A  Search  in  the  Fimction  Space  by  Directed  Genetic 
Programming,  D.  E.  Makarov  and  H.  Metiu,  J.  Chem.  Phys.  108,  590  (1998) 

5.  Electrons  Solvated  in  Zeolites,  N.  P.  Blake  and  H.  Metiu,  Proc.  Australian  Summer  School  on 
Condensed  Matter  Physics  (January  1997),  M.  K.  Das,  editor.  World  Scientific,  Singapore  (1998) 

6.  Migration  of  Hydrogen  on  a  Solid  Surface:  The  Physics  of  the  Process  and  the  Methodology, 
H.  Metiu,  Proc.  Fermi  School  on  the  Dynamics  of  Rare  Events  (Lerici,  July  1997) 

7.  The  Rate  of  Photon  Absorption,  H.  Metiu,  Proc.  Fermi  School  on  the  Dynamics  of  Rare 
Events  (Lerici,  July  1997) 
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Computation  of  a  Rate  Constant  and  Their  Use  in  Improving  and  Testing  the  Quality  of 
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1.  Other  funding 

There  is  no  relationship  between  the  work  supported  by  this  grant  and  the  following: 


National  Science  Foundation  Chemistry  Division  “Theoretical  Studies  in  Quantum  Dynamics” 
8/1/94-7/31/98,  $365,000  total 

Air  Force  Office  of  Scientific  Research“Computer  Simulation  of  Growth  Phenomena” 

1 1/1/94-10/31/97,  $303,500  total 

National  Science  Foundation  “Center  for  Quantized  Electronic  Structures  (QUEST)”,  1993-98 
(one  of  17  co-PIs) 
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a.  Principal  Investigator:  Horia  Metiu 

b.  Current  telephone  number:  (805)  893-2256 

c.  Cognizant  ONR  Program  Officer:  John  C.  Pazik 

d.  Program  objective 

We  use  theoretical  methods  to  help  Prof.  Galen  Stucky's  group  in  their  search  for  better 
thermoelectric  materials.  We  examine  a  class  of  compounds  based  on  Si,  Ge,  or  Sn  which  have  a 
clathrate  structures.  There  are  about  100  such  compounds,  each  allowing  a  fair  number  of 
chemical  variations  (doping,  slightly  non-stoechiometric  composition,  etc).  They  are  difficult  to 
synthesize  and  characterize  structurally  and  the  measurements  of  their  transport  properties  are 
very  tedious.  For  this  reason  it  would  be  extremely  useful  to  provide  theoretical  guidance  and 
develop  methodology  that  will  allow  us  to  determine,  before  synthesis,  which  among  these 
compoimds  is  likely  to  have  the  best  thermoelectric  properties.  We  are  providing  such  guidance 
by  performing  density  functional  calculations  of  the  band  structure  of  these  compounds  and 
using  it  to  estimate  the  Seebeck  coefficient  and  the  electrical  conductivity.  These  calculations, 
though  difficult,  are  more  rapid  than  experimental  synthesis  and  characterization.  Our  work  will 
allow  the  experimentalists  to  focus  on  the  more  hopeful  candidates.  We  are  also  proposing  the 
synthesis  of  new  clathrates  based  on  GaAs,  AlAs,  In  As,  CdTe  and  other  tetrahedrally  bonded 
compounds.  By  performing  total  energy  calculations  we  can  determine  which  among  these 
compounds  are  stable,  as  compared  to  other  solids  that  can  be  made  from  the  same  atoms,  and 
therefore  have  a  good  chance  of  being  prepared  in  the  laboratory. 

e.  Significant  results  during  last  year 

Our  main  effort,  in  collaboration  with  Prof.  Galen  Stucky,  has  been  directed  towards  determining 
whether  zeolites  containing  one  solvated  electron  in  each  cage  are  good  thermoelectric  materials. 
For  this  we  had  to  understand  the  electronic  structure  of  the  material.  Calculations  that  explained 
most  of  the  measured  quantities  indicated  that  the  material  is  a  metal  with  low  electron  density. 
However,  disagreement  with  the  low-temperature  magnetic  susceptibility  measurements  forced 
us  to  substantially  refine  the  model.  By  using  density  functional  calculations,  which  include 
generalized  gradient  and  self-interaction  corrections,  we  have  shown  that  the  system  is  an 
antiferromagnetic  Mott  insulator.!  This  means  that,  in  spite  of  the  fact  that  the  Seebeck 
coefficient  and  thermal  conductivity  are  favorable,  this  class  of  new  substances  will  not  contain 
good  thermoelectric  materials  since  they  vvill  not  have  sufficient  electrical  conductivity.!. 2  In  this 
work  we  have  reached  a  very  important  conclusion:  it  is  not  possible  to  obtain  even  a  qualitative 


4 


understanding  of  the  electronic  properties  of  a  narrow-band  solid  without  making  the  self¬ 
interaction  correction,  which  is  usually  neglected  in  solid  state  calculations.  Another  important 
achievement  is  the  development  of  a  new  method,  directed  genetic  programming,^  to  find 
potential  energy  functions  that  fit  the  energies  calculated  by  ab  initio  methods.  This  difficult, 
tedious,  and  important  problem  can  now  be  given  to  computers. 

I  was  invited  to  lecture  on  electrons  in  zeolites  at  two  summer  schools  (in  Australia^  and 
Europe^’^)  and  this  gave  me  a  chance  to  review  some  of  our  ONR-sponsored  work  for  a  broader 
audience  of  physicists  and  chemists. 

The  bulk  of  our  work  this  year  has  been  developing  software  for  calculating  the  electronic 
and  transport  properties  of  semiconductor  clathrates  (see  item  f). 

f.  Brief  summary  of  plans  for  next  year’s  work 

We  are  now  concentrating  on  class  of  clathrate  compounds  that  promise  to  be  excellent 
thermoeleetrics.  There  are  four  types  of  interesting  clathrates,  described  by  the  formulae 

^8^48’  A8Y8T38,  and  BgYi6T38,  where  T  are  tetrahedrally  coordinated  atoms  (Si,  Ge, 
Sn);  A  are  alkali  atoms;  B  are  Sr,  Ba,  Ca;  and  Y  are  Al,  Ga  or  In.  So  far  the  transport  properties 
of  only  one  compound  of  each  type  have  been  studied  experimentally  and  they  all  have 
promising  thermoelectric  properties.  One  of  them,  Sr8Gai5Ge38,  has  a  figure  of  merit  half  that  of 
the  best  commercial  compound.  There  are  about  100  clathrates  and  their  transport  properties  can 
be  changed  by  chemical  modifications.  It  is  extremely  likely  that  some  among  them  are  better 
thermoelectrics  than  the  ones  already  studied.  We  plan  to  use  density  functional  band  structure 
calculations  to  aid  the  experimentalists  find  which  clathrates  have  a  chance  of  being  the  best 
thermoelectrics.  In  these  calculations  we  determine  the  band  structure,  and  estimate  the  Seebeck 
coefficient  and  the  electrical  conductivity.  We  have  spent  a  lot  of  time  developing  the  software 
needed  for  these  calculations  and  have  just  started  performing  them.  We  have  found  that  S^g  is  a 
semiconduetor  with  a  larger  band  gap  than  Si.  NagSi4g  is  a  metal  formed  because  the  alkali  atoms 
donate  electrons  in  the  conduction  band  of  Si4g.  The  compound  NagAl23Si23  is  a  metal  and  it  is 
likely  to  have  a  low  Seebeck  coefficient.  Such  calculations,  while  extremely  laborious,  are  faster 
than  the  synthesis  of  compounds  and  measurements  of  transport  properties  and  will  serve  as  a 
guide  to  pick  the  promising  compounds  for  synthesis. 

g.  List  of  names  of  graduate  students  and  post-doctoral(s)  currently  working  on  the 
project:  Dr.  Paul  Weakliem,  Dr.  Nicholas  Blake,  Dr.  Dmitrii  Makarov 


PARTIU 

Figure  1  (5-part).  We  show  the  structure  of  a  clathrate  having  48  frame  atoms  in  the  unit  cell. 
The  frame  atoms  ( e.g.  Si  or  Ge)  make  cages  which  are  filled  with  alkali  or  alkali  earth  metals. 


Figure  2  (on  page  2).  The  sodalite  cage  consisting  of  -Si-O-Al-O-Si-O-Al-O-Si-  chains.  The  solid 
are  made  by  joining  the  cages  at  their  hexagonal  faces.  There  are  four  ions  and  an  electron  inside 
each  cage. 


Figure  3.  There  is  no  figure  3. 


6 


s  a 

•2  .2 

o  o 
C3  c^ 

M  u 

«S  B 

rt  ^  *3 

W  •r-^  Cj. 

O  c/5  23 

•JS  d  (D 

.ts  w  t/i 

««  'O  <-! 


e«  ^  S 

O  .2 

t3 

«  c  3 
&  U  S 
S  ^13 

HH  •  O 


r^  O 

w  3 
§  o 
^  i:3 
2  ^ 
ia  CO 
^  <L>  3 

3  sa  ‘c 

6  :S  3 

§  3  S 

•-u  +3  0) 

■K  CO  CO 

sa  *-<  w 

<l>  (1> 

^  -a  ’3 

O  S=l  c+_, 

0^3  0 


bO 

•2  d 

ti  3  o 

^  -8 

I  8 

4->  ^ 

•  r-H  0  1-0 

3  d  o 

(D  ^  O 
N  bO 

TS  ^  <L) 

53 

-Q  c3  S 

;>.  5^  C> 
►O  Vh  -M 


^  S<  jl3  CO 

‘S  S  ^  8 

1 — I  z3  5 

^  g  d  « 

•  .S  (D  'd 


(U 

2 

^  8) 

4h  .. 

-2 

fi  g  .3 
o  ^  U 

X  ro  C-/ 

&  <u 

S  CO  13 

O  So 

§  >  g 

•  :zi  ’B 


O  rt' 

Q 1 

O  o 

CfH 

CO  o 
<D  ^ 

•rH  -K 

ti:  i 

(D  -a 
Ph^ 

o  .a" 
PM  s 


c^-l 

d 

O 

O  •  r! 

3 

CD  5-1 

'o  <l-> 

«  -S' 
^  ‘fc! 

^  (D 
(D  Vh 

"  c 

(u  S 

cd  >^ 
(U  ^ 

^  <u 

c  o 

o 

o  2 

B  5i 

Sf— < 

<D 

^  S 

HH  ^ 

^<§ 
>3  CO 

S  'S 

pG  0 

a>  cd 

H  S 


bX) 

Ch  'd 
O  (D 
g  Ph 

s  -§ 

o  ^ 
o 

<u  ^ 

^  CO 

O  CO 


s? 

•  f-h  0^ 

■s  ^ 

s  a 


^  C  ti 
<D  G  cS 

S  ^ 

'S  G  O 

c2  "S  2 

C4H  O  d 

»G  CO  3 
O  d  G 
CO  o 

G  CO  Q 
D  d 

^  13  B 

O  £H  P 

.5  CO 
lo  CO  -d 

"33 

>  G  G 
T3  G  CO 
O  O  (D 

4)  .d 

0  1— H 

*p  S  P 

^  o  d 
O  N 


d  -d 

^  • 


C+H  O 

3 

CO  t! 

-S  ^ 

'>  ^ 

^  .2 

o  ^ 

d  P 

d  o 

B  « 

-d  ^ 

a  o 
G  -d 
.2  P 

*13  c3 

^  la 

-H  0^ 

cB 


cj  d 

B  s 


(D 

Pi 

4I| 

-i-> 

<U 

<u 

'S 

Td 

o 

(U 

B 

3 

13 

§ 

s 

'Trt 

o 

§ 

G 

o 

•  t-H 

4-^ 

o 

"Td 

•  »-H 

o 

o 

C 

U 

(U 

o 

•  T»^ 

CO 

o 

> 

C+-H 

•  t-H 

CO 

a> 

t3 

o 

CD 

X) 

'"S 

nd 

c 

(U 

3 

<D 

(U 

a 

o 

nd 

3 

00 

o 

CO 

O 

•  r- < 

O  N 


a  M 
a  G 


^  d 

•§  §^ 

2  bO  d 

^  2  G 


d  d 
'G  rH  « 

d  S  ^ 

Cm  O  '5 


C3  d 

G  .s 
•2  9 

o  ro 
G  ^ 


d  CO 

w  d 
G  O 
d  d 
d  3 

g  o 

^  d 

•  o 


^  d 
d  d  o 
to  'G 

'2  B  3 

^  ■§'8 

bO'rG  ^ 

CO  (L> 
Q  ri 

Cj  CO 

S  ?-i 
2  ^  .2 
O  ^ 

d  3  & 
^39 


CO  CO 

<D  (U 

•  1-H  4-» 

tlj 

O  cd  0 
^3  G 

p  3 
B  .BB 

-2  2  d 

I  s  « 

^|:s 

C4H  o 
•  0-2 


3  t: 

3  d 
d  Gh 
d  o 
>  G 
d  Cl, 


1/3 

a 

'M  ’S 

«4H  5 

^  ■ 
^  A 

«  o 

a  -3 

i4 

^  /A 
&x  O 
a  « 

f  » 

o  ^ 

^  g 


IZl  O) 

T.  ^ 

•  pH  W 

HM  _ _ 

as  •£: 

t3  o 
o  « 

C/3  W 


1/3 

^  s 

3  is 

S  U 
‘S  i> 

fl 

G  S 

O  ^ 
WO  S2 

««  5 


-C  03 

pS  ^ 
ts 

'Sod 

p  ^ 


o  <1> 
•!=4  WD 

oj  Cd 

£5  ^ 

Cd  jCS 
C/D  U 

a.  s 
s  ^ 

^  fi 

H  S 


la  c^ 
V3  O 

Si  I 
pP  p 

fe  I 
S-3 

ns  'O 
S  © 

N  !P 
P3  P 
P  « 

©  1/3 

p2 

2  '3 

'P  *2 

</)  4^ 

•  PH  ^ 

p  2 

®  S 

t3  ^ 

Si  Si 

p 

<u  .2 

pp  pp 

H  H 


Vi  Vi 

2  .Si 

tS  ;h 

*©  ^ 
O  4) 

O  § 

^  g 

•M  c; 

OJ  pP 

-p  ^ 

'M  ns 

5  © 

a  w) 

©  Si 
^  p 

©  Vi 

TS  TJ 
^  P 
>  P 

P  © 
pP  ft 

©  s 

^  8 


p 

pp  p 
p 

^  a 

’■§  i 

u  h 
ft  tj 

Si 

^  ft 

P  Vi 
p 

S  O 

T3  '5 

s  & 

I  8 

^  -S 

«M  P 
Od 

pT.Sf 

©  ^ 

H-»  CJ 
C/D  ^ 

cu  u 

rj^  ^ 
© 

-p  © 

TS  Vi 

■P  P 
P  ^ 

Vi  H-i 

P  . 
.2  « 
p  'S 

p  S 

-2  ^ 

^  •- 

©  1/3 

P^  P 
P  © 
P  'P 

®  Si- 

^  o 
2 

®  ea 

«M  t« 

.f  i 

s  ^ 
o  ^ 


p 

•  5rt 

©  © 


'2  " 

p  ^• 
^  2 

2  © 


2  c 

t3  p 


©  »-' 

2 


^  -I 
£  © 


u 

-s  a 

P  t/3 

s  § 

^  *-c 

^  p 


P  -PH 

2 


vrj 

^  CQ 

a 

H-d  T5 
^  2 
u  u 

©  P 
© 

2  © 
©  © 
-*-»  ^ 
T3  wo 
P  © 
p©  P 

^  Vi 

p  p 

2  ® 
'p  tp 
U 

^  O) 

C  ^ 

®  fe 

•r^  W 


P  *P 
-  H 

-li  c 
2  .2 
§  tl 

1  ^ 

p  © 

^  P 

^  © 
•P  *X3 

2  © 

I  2 

^  <u 

•S  p 

S  •fn 

I 

s?  a=i 


•  c^ 

p  P 

Vi  p 

c«  5 

8  2 

^  wo 

c  ^ 

C  'tri 

"  s 

•ti  wo 

Vi  P 


P  ^ 

>-'  4; 

p  SP 
®  ti 
H  P 

tS  «« 

a  ti 

©  © 

S  8 

c  ^ 

.2  p 

p  p 
N  "P 
PS  .p 
p  ® 

«  I 

©  © 

^  S 

p  ® 

-p  ."p 

^  .HH 

«M  P 

©  5 

P3  2 

©  •P 
•p  .P 
©  ^ 

S  p 
a 

Si  ^ 

p  p 
p  2 
^  «2 
43  © 

•N 

pg  ^ 

PS 

®  ^ 
«/3  ^ 


(Jnce  the  correction  is  made  it  changes  qualitatively  the  nature  ot  the  material:  insteai 
of  being  a  metal  it  is  a  Mott  insulator.  This  means  that  no  doped  zeolite  will  have 
enough  electrical  conductivity  to  be  a  good  thermoelectric.. 


A  theoretical  study  of  clathrates 

The  thermoelectric  properties  of  only  four,  out  of  about  a  hundred  clathrates,  have  been 
measured  experimentally.  They  are  all  good  thermoelectrics  and  one  (Sr8Gai6Ge3o)  is 
close  to  the  best  commercially  used  compound. 
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in  their  properties:  Si4  ^  is  a  wide  band  gap  semiconductor.  Nag  Si45  is  a  doped  Si46  with 
metallic  character,  Na8Al23Si23  is  a  metal  and  Kg  Ge4g  is  a  semiconductor  and  a  good 
thermoelectric  material.  This  variability  gives  us  hope  that  we  can  find,  among  a 
hundred  of  such  compounds  and  all  their  modification,  a  superior  thermoelectric. 


